Quasi-polynomial mixing of the 2D stochastic Ising model with "plus" boundary up to criticality

We considerably improve upon the recent result of [37] on the mixing time of Glauber dynamics for the 2D Ising model in a box of side L at low temperature and with random boundary conditions whose distribution P stochastically dominates the extremal plus phase. An important special case is when P is concentrated on the homogeneous all-plus configuration, where the mixing time T-MIX is conjectured to be polynomial in L. In [37] it was shown that for a large enough inverse temperature beta and any epsilon > 0 there exists c = c(beta, epsilon) such that lim(L ->infinity) P(T-MIX >= exp(cL(epsilon))) = 0. In particular, for the all-plus boundary conditions and beta large enough, T-MIX <= exp(cL(epsilon)). Here we show that the same conclusions hold for all beta larger than the critical value beta(c) and with exp(cL(epsilon)) replaced by L-c log L (i.e. quasi-polynomial mixing). The key point is a modification of the inductive scheme of [37] together with refined equilibrium estimates that hold up to criticality, obtained via duality and random-line representation tools for the Ising model. In particular, we establish new precise bounds on the law of Peierls contours which complement the Brownian bridge picture established e.g. in [20, 22, 23]